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1. Introduction 

Consider classical Yang- Mills theory. If the initial conditions contain field energy in some 
region of space with vacuum around, we know that after a while the energy will be in an 
expanding nearly-spherical shell of outward moving waves. But if the initial conditions are 
statistically spatially homogeneous, the fields will re-interact indefinitely. Now classical 
Yang- Mills theory, like all nonlinear classical field theories, has no equilibrium; there is in- 
finite phase space in the deep ultraviolet. So the energy should cascade into the ultraviolet, 
spreading with ever smaller amplitude over ever more field modes. But exactly how does 
this cascade actually occur? 

Recently there has been renewed interest in this question, because classical Yang-Mills 
theory should describe the behavior of weakly-coupled, quantum Yang-Mills theory when 
the mean occupancy is high, f{k) ^ 1 for the most important in a system. The initial 
conditions after a heavy ion collision, in the limit of large and high energy nuclei, are 
expected to involve high-occupancy, effectively classical fields |l|]. It is not clear to us how 
realistic the classical field description is phenomenologically; also, the initial conditions for 
a heavy ion collision are spatially inhomogeneous and anisotropic. Nevertheless, we will 
consider the homogeneous and isotropic case as an interesting warmup problem and lesson 



- 1 - 



in the behavior of classical Yang-Mills theory. At minimum we believe that it is interesting 
to study the behavior of classical Yang-Mills theory in its own right, especially to establish 
its differences from scalar field theory. 

Quite recently two papers have considered essentially the problem we raise; Yang-Mills 
theory at weak coupling and with initial conditions with high enough typical occupancy that 
the behavior is that of classical fields (at least for some period of time). Two of us argued 
m that the cascade into the ultraviolet would proceed with the dominant momentum scale 
Pmax (the scale where most of the energy density resides) growing with time as Pmax oct^. 

4 

The occupancy at this scale would behave as f{pma.x,i) oc t~T , and the occupancy for 
P ^ Pmax would scalc as f{p) oc p^^. (The same scaling for Pmax was found for scalar field 
theory much earlier by Micha and Tkachev Q.) On the same day, a paper by Blaizot, 
Gelis, Liao, McLerran and Venugopalan [Q] came to similar conclusions, except that they 
argued that there could in addition be the formation of a condensate of infrared excitations, 
carrying most of the particle number, but a minority of the energy density, in the system. 
A recent numerical study of the problem by Berges, Sexty and Schlichting HI did not 
try to determine the scaling of momentum or occupancy with time, but investigated the 
infrared tail of the spectrum. Contrary to the arguments of they found f{p) oc 

with a = 3/2 at early times and transforming to a = 4/3 at late times (an exponent first 
proposed in Reference @). They found no direct evidence for a condensate, but argued 
that the a = 3/2 scaling might be indicative of a condensate's physical effects. 

We will re-examine this problem, using lattice methods. We begin by reviewing the 
arguments for Pmax ^ 

1 7 and f (xt 7 scaling, laid out in ||2|, Q . Then we study the evolution 
directly on the lattice by initializing large volume lattice systems with random, infrared- 
dominated and nonlinear-amplitude initial conditions. We study their time evolution over 
long times, and with detailed control of lattice and initial condition effects. We verify that 

1 —4 

the 

Pmax OC and f (X t T scaling is established very fast, but the infrared behavior with 

_ 4 

the p 3 scaling emerges somewhat slower; and the interpretation of the infrared region 
is complicated by the physics of screening. Once the physics of screening is taken into 
account, we see no evidence for infrared scaling exponents other than —4/3. We do not see 

-4 

infrared occupancy in excess of that predicted by a p 3 exponent, but we nevertheless study 
what would be the properties and fate of electric (plasmon) and of magnetic condensates 
in the infrared; in each case we find that they would be short-lived. 



2. Scaling in Classical Gauge Theory 

We will start by explaining the relation between classical Yang-Mills theory and the usual 
quantum theory at weak coupling, and we will introduce the scales in the problem. Then we 
see how the scaling behavior naturally arises by considering the structure of the scattering 
terms in the Boltzmann equations. This section reviews known arguments, see for instance 
Refs [0, 0]; readers familiar with its contents may want to skip this section, except for 



Eq. ( |2.2| ) and Eq. (2.12), which define our scale Q and momentum p„ 
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2.1 Classical Yang-Mills theory and scales 

Without h, length and energy scales are distinct, and the gauge coupling is dimensionful. 
The easiest way to think about the classical theory, for someone familiar with the quantum 
one, is to assume all occupancies are of order f(p) ~ 1/(7^, so f{p) is a number of order 
unity. Then one takes small (to zero) holding f{p) fixed. This immediately means we 
can ignore fermions; since we assume no scalar matter, we will work with pure-glue QCD. 
Occupancies are naturally /(p) ~ l/s'^- If Q is a scale characterizing the wave- vectors of 
initial fluctuations, then the energy density is naturally e ~ Q'^ /g'^\ it is better to interpret 
Q,p as wave numbers or inverse lengths than as energy or momentum scales, (l/y^ plays 
the role, dimensionally, of h, turning an inverse length Q into an energy scale Q/g'^ ) When 
9^ f{p) is small, the theory behaves nearly linearly; when g'^ f{p) ^ 1 the behavior is highly 
nonlinear, and the notions of wave- vector and occupancy become highly gauge dependent. 
Here "small" does not mean oc g"^ , since we take g"^ small first. Rather, the behavior is 
perturbative when there is some other expansion parameter A, with f{p) ~ A/^^, and 
A <C 1 but g^. We will see that the inverse system age - actually t~7 - plays the role 
of A. 

We will consider extensive systems with a mean energy density e ^ Q'^/g'^. To make 
quantitative statements easier, we will define the scale Q in terms of e as follows. When 
Yang-Mills theory is nearly linear, it makes sense to fix the gauge and describe it in terms of 
quasiparticle excitations with occupancy /(p), in terms of which the energy density would 
be^ 

s'jVcE ^ / k> [g'NJlk)] dk . (2.1) 

where the sum is a sum over spin and color states, J2sc ^ ~ ^(-/V^ — 1). With this in 
mind, we simply define the scale Q as 

,^.,.Mizl)g. g.(^y (2,2, 
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SO that, in the perturbative regime, we have 

~ j g'^Ncfik) dk . (2.3) 
We will generally use dimensionless time Qt and dimensionless wave numbers p/Q- 
2.2 Scaling from the Boltzmann equation 

Now consider the Boltzmann equation describing the time evolution of the occupancy f{p). 
According to Arnold, Moore, and Yaffe 0, it is of generic form 

df{p,t) 



dt 



-C2^2[f{p)]-Ci^2[fip)] 



^Here we assume we are considering SU(A'c) gauge theory; for a general gauge group, replace (A'^^ — 1) with 
the dimension of the adjoint representation dj^ and A'c with the second Casimir of the adjoint representation 
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C202[/(P)] = ^ [ \Mpk;p'k'\^ (2T:)H\p+k-p'-k 
Jk;p'k' 



f{p)f{m+f{p')][i+f{k')] - [i+f{p)][i+f{k)]f{p')f{k')) (2.4) 

' — r c^^fc 

k J (27r)^2fc" ■ 



(we discuss Ci^2 in a moment). Here the phase space integrals are = / (2n)Hk^ • ■'^^^ 



pure-glue QCD we have 

|AT..;,'.'P = 4iVcV(3-f -^-^) (2.5) 

with s,t,u the usual Mandelstam variables, and with the su/t^ and st/u^ matrix elements 
IR regulated as described in Ref ||7|. We are interested in the behavior for / ~ 1/(7^ 1, 



which allows a little simplification. The terms in Eq. ( p.4D proportional to cancel 
between "loss" and "gain" terms (the first and second terms in the large parenthesis in the 
last line); but the terms of order do not; and since ^ we can drop the terms 
and write 



[mfim^+fip'm+fik')] - [i+f{pm+m]f{p')f{k')) 
- f{p)f{k)f{p')f{k')[r\p') + r\k') - r\p) - r\k)) for / » i . (2.6) 

For generic values of / and for large momentum transfers^ so s ~ — t ~ —u, there are 
no cancellations in the round bracket in the last line, and this term is of order Next, 



Eq. (2.5) shows that \M\ is of order g . Therefore both sides of Eq. (2.4) are of order 
g^"^ . Since |A^P is of order N^g'^, it is natural to assume f{p) oc l/{g'^Nc). If we introduce 
/ = g'^Ncf, then all factors of g'^Nc cancel when we express the Boltzmann equations 
in terms of /. An explicit factor of g'^Nc also cancels the factor associated with f{p) in 
determining the screening mass, needed in the IR regulation of A^^. Since the leading 
order expressions for Ci^2 and C2<->.2 do not contain subleading-in-l/A'^c dependence, there 
is also no dependence on Nc left, which means that a study using the group SU(2) should 
return the same scaling behavior as any other SU(A'c). 

With g^ taken care of, we now examine the expected scaling properties with time. 
It is natural to assume that, a time t after some initial conditions are established, and 
considering a typical momentum in the range which dominates the system energy density, 
the time derivative of the occupancy is of order^ 

df{p,t) f 

—dT- ^ 1 ■ ^'-^^ 

Since the matrix element is dimensionless, in terms of the (time dependent) typical mo- 
mentum scale Pmaxit), the momentum and /-scaling of the two sides of the Boltzmann 



^When p' ~ p so — t ^ s, diverges as s^/t^ but there are compensating cancellations in Eq. (2.6), 
rendering the coUision term at worst logarithmically divergent, cut off by screening efltects. 

^If df /dt 3> f /t then the occupancy quickly adjusts to be something different; and df /dt <ti f /t can only 
occur at early times for low-occupancy initial conditions. If we are interested in late times, the absence of 
a thermal ensemble ensures that the occupancies will eventually show evolution with df /dt ~ f /t. 
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equation are 

(2.8) 



t 

This has a simple interpretation; the scattering rate T for a typical particle should be 
Tt ~ 1. The scattering rate is 0{g^praaxf^) = C'(Pmax/^); the factor Pmax is on dimensional 
grounds, g'^ is because there are two vertices involved in a 2^-7-2 scattering process, and 
is the occupancy of the scattering target and the stimulation factor for the final-state 
scattering target. (The stimulation factor for the outgoing particle under consideration 
cancels between gain and loss terms.) The Ci^2 term has the same g'^ and / scaling.^ 
Furthermore, energy conservation implies that 



g'N.e -Q= J r /(P) dp ~ pLx(i)/(Pmax, t) pLx(i)/(Pma., t) ~ (2.9) 

and in particular this combination is time independent. Solving these two relations, 
Eq. ( |2.8| ) and Eq. (|2.9|), for / and Pmax in terms of Q and t, we find 

g(Qt)^ , (2.10) 

/(Pmax,t) ~ (Qt)"^ (2.11) 

the scaling behaviors already found in Refs [§, (A of Ref. @] refers to the same parametric 
scale as Pmax)- 

Since we will use Pmax 

repeatedly in the following sections, it behooves us to define it 
more precisely. We will define it as 

2 _ {TrjDxB)^) ^ fk'f{k)dk 

Pmax - 1 + ^2)) " J ^ f{k) dk ' ^ ' ' 

where the angular brackets indicate volume (and/or ensemble) averaging. The definition 
here is convenient because it is gauge invariant and easy to evaluate on the lattice; the 
second expression is the small-amplitude or quasiparticle behavior of the first expression. 

The arguments above suggest the same scaling for scalar fields should also occur Q. 
But for scalars there is an extra complication. The total particle number density is n ~ 
/ fip)p^dp, which according to Eq. (|2.10| ) and Eq. ( 2.11 ) should scale as n oc t~ . But 



C202 does not change particle number, so no scaling solution can actually solve Eq. (^ 
if only C2<->2 is present. In a scalar theory this means that the extra particle number must 
somehow crowd into the infrared, where it either forms a condensate or is destroyed by 
nonlinear processes, which occur faster at the very high occupancies which will occur in 
the infrared (for a recent study see Ref. |8|). For us, Ci<^2 is of the same order as C2^2 5 so 
particle number is not conserved and there need not be a condensate or enhanced infrared 
occupancy. Whether or not such infrared enhancements occur in practice requires a more 
detailed solution of the problem, which we turn to next. 



*Very briefly, the l-f->-2 process involves a smaU-angle scattering, with a sphtting. The rate for small-angle 
scattering is g^pf^ like €2^2, times a soft enhancement /m^, which is the ratio of the typical momentum 
to the screening scale rn^ ~ g^P^f- And the splitting introduces an additional factor of g^f, with / the 
stimulation associated with the extra particle. Putting it together, the rate is g'^pf'^ip^ /{g'^p^ f)){g'^ f) ~ 
g'^pf^ , the same as the rate of large-angle 2-f->2 processes. 
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3. Lattice treatment 



We will solve the Boltzmann equations for this system in a future publication. Here we 
will directly solve classical Yang-Mills theory fully nonperturbatively, on the lattice. Our 
treatment will emphasize finding the scaling solution and identifying (transient) corrections 
to scaling. We will also make an effort to investigate the possibility of infrared condensates. 
There are large algorithmic advantages to considering only SU(2) gauge theory, so this 
study will strictly work within SU(2). 

The technology for studying real-time, classical fields on the lattice is decades old; see 
for instance |p. The main issues are, what initial conditions should be studied and how 
much do results depend on the choice, what physical measurables should be used, and how 
does one monitor finite volume and finite lattice spacing effects? 

|— I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — I I — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — I 



Figure 1: Evolution of the characteristic momentum scale p^ax * function of time t, each scaled 
by Q to be dimensionless. At the left, time is based on the starting time of each simulation; at the 
right, the starting times are shifted so the late-time behavior is in better correspondence. In the 
righthand figure, the red dotted line is a fit based on strict p^^x {Q'^)^ scaling behavior. 

The simplest gauge invariant measurable is p^ax defined in Eq. ( ^.12|) . So to examine 
the approach to scaling we will first look at results for this quantity. To do so we created 
several different initial conditions, with either infrared electric fields, magnetic fields, or 
both present, and with different ranges of momenta occupied. To keep lattice systematics 
common, they all have approximately the same value of Qa, namely Qa ~ 0.42. The 
evolution behavior of p^^^/Q^ ^ ^ function of Qt is shown for 6 initial conditions in Figure 
|l[ Each curve shows an initial transient, whose length depends on the initial condition; but 
after some time they all approach a common scaling behavior. In the lefthand figure the 
time is based on the starting time of the simulation; but depending on the initial condition, 
there may be a delay before the dynamics starts to track towards the scaling solution (or 
it may reach the scaling solution early, if the initial condition is already quite similar). 
Therefore, to understand whether the evolution really approaches a scaling solution, it 
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Time Qt 
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makes more sense to shift each initial time so that the curves are more similar at late time. 
We have done so in the righthand figure, which shows that the different initial conditions 
very accurately fall onto the same scaling behavior. The red dotted line is a 1-parameter 

2 

fit, 

Pmax(^) — (^(Qt)"^ , which shows that the curves are obeying the expected time scaling. 
Unfortunately, by Qt = 450 the scale Pmax has already reached Pmax« = !> which means 
that lattice spacing errors start to occur; at Pmax« = 1 we expect of order 10% errors due 
to higher-dimension operators in the lattice action. So precision quantitative results would 
demand a finer lattice (smaller Qa). 

1 n — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 2: Ratio of electric to magnetic field energy for the same evolutions as in Figm'c ^ Initial 
transients rapidly decay onto a tracking solution with slightly more electric than magnetic energy. 

We can also look at the gauge invariant ratio of electric to magnetic energy densities 
(Tr E^)/{Tr B^), shown in Fi gure 0. A quasiparticle with momentum large compared 
to the plasma frequency p ^ Wpi should have equal time-averaged electric and magnetic 
energy. Plasmons with p < Upi would display more electric than magnetic energy, whereas 
magnetic fields in this regime (the Landau cut) would have more magnetic energy. So a 
condensate of plasmon quasiparticles would be expected to manifest as an excess of electric 
field energy leading to a ratio larger than 1, while a condensate of magnetic field would 
give a ratio smaller than 1. In fact, after a short and initial-condition dependent transient, 
the ratio settles down to be close to but slightly above 1, decaying towards 1 with time.^ 
This suggests against the existence of large condensates; but we will examine the evidence 
in more detail in the next section. 



3.1 Gauge fixed observables 

To learn more about the evolving configurations on the lattice, it is useful to examine gauge- 
fixed observables. Coulomb gauge is the choice of gauge which minimizes It is 
therefore perhaps the most sensible gauge to use for studying the distribution of excitations 
at one moment in time. And the technology for fixing to Coulomb gauge on the lattice 

^Specifically, in the figure, the blue curve has initial conditions with purely electric energy; the yellow 
curve has purely magnetic energy. The other curves have a mix of electric and magnetic but in different 
momentum ranges. 
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is very well known, amounting to the fixing of Landau gauge for the 3-D configuration at 
an instant [^]. Since the gauge fixing procedure emphasizes minimizing A^(x) at every 
time at the expense of minimizing A^, unequal time correlators in this gauge tend to have 
very short autocorrelation. So we will not attempt to untangle unequal time correlations 
using gauge fixed methods. 

Conventionally, one uses the perturbative relation between occupancy and field ampli- 
tude, for transverse quasiparticle excitations, 

I d'x e^^-{Al{x)Aim = ^^^^^/(P) , (3.1) 
I d'x e^P-''{Elix)EiiO)) = (5abVl}{p) \p\)f{p) (3.2) 

(with Vl^ = 6^^ — p^fP the transverse projector), to estimate the occupancies as 

fA{p) = ^^^IpI / e^^-(AUx)Al(0))eoui , (3.3) 

fE{p) = 2[Nf-l)\p\ I '^''^ ^'^■'(^a(^)i?^(0))coui . (3.4) 

In this way we obtain two estimates of the occupancies, one from gauge fields and one 
from electric fields. To illustrate the utility, and the danger, of this approach, we use it to 
examine the occupancies of a nonperturbative lattice system in equilibrium.^ In equilibrium 
one usually expects f{p) = T/p, which is small provided p <^ g^N^T. To ensure that the 
occupancy is small for typical lattice momenta p ^ ii/a, we consider a lattice with spacing 
a = OA/{g'^ NcT) (which is /3 = 20 in the conventional lattice notation). The results for 
the occupancies Ja and /e are shown in Figure |^. 

The figure shows /a and /e as the black and green curves. But whereas the ^-field is 
automatically transverse because of our gauge condition, the electric field is not. Gauss's 
law states D ■ E = 0, not V ■ E = 0, and below the Debye scale the two become strongly 
inequivalent. Therefore we have separately plotted the transverse electric estimator 

JeM = 2{Nf-l)\p\ I ^'"^ ^^^■"(^?(^)^'(0))coui (3.5) 
and the longitudinal electric estimator 

^^'^^^ = (i^-l)lpl / e^^-"(^f (^)^i(0))coui . (3.6) 
At the level of hard thermal loops these have simple expressions (see Appendix ^ 

m? T T 
Ie, (p) = 2 ° 2 - ' fEt iP) = -- (3.7) 

TTlp + p"^ p p 



®The lattice system has an equilibrium because the available momenta are restricted to a Brillouin zone, 
so the phase space is finite. From the point of view of continuum classical field theory, this equilibrium 
state is a complete lattice artifact. 



-8- 




Dotted:model . 



2 4 6 

Momentum k/(g2NJ) 

Figure 3: Equilibrium occupancies in a 256'^ SU(2) lattice with a = 0.4/ {g2NcT), scaled by k/T so 
at lowest order they should be 1. The black curve is /a estimated using the ^-field correlators; green 
is Je estimated using i?- field correlators, while red and blue are the estimates using the transverse 
and longitudinal components of the i?-field. The most faithful estimator is the transverse E-field 
component. 



The estimator /e is /e = fEt~^\fEi, which therefore systematically rises at small momenta 
p ^ mi). Physically, the reason /e rises is that, for p < niD, there are three polarization 



states contributing to the electric field correlator, but the estimator Eq. (3.4) assumes 
that there are only two. Therefore, below the scale mj), the estimator /e is incorrectly 
normalized and systematically high. However the estimator fEt is based on a correct 
counting of degrees of freedom and should be used instead. 

The figure also shows that the >l-field correlator does not reproduce the expected 
behavior fAik) = T/k. As one goes down in momentum scale k, it first rises above the 
leading-order estimate, peaks, and falls to zero at small momentum. The rise is due to 
perturbative corrections to the AA correlator, which are known; at next-to-leading order, 

p — 6f-\p\ + c>y9 T^) 

Note that the first correction is suppressed by 1/p, not by 1/p^. The dotted "model" curve 
shown in Figure ^ uses the above correlator with 0{g'^T'^) replaced with g^N^T"^ /36, which 
seems to match the infrared behavior. The "model" curves for the Je curves are based on 



Eq. ( |3.7| ) and the lowest-order lattice value for m^, [11| 



m\,c? = — '!— = -^ , S = 3.1759.... (3.9) 
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Because /a rises above the leading-order expectation, using / 0fA{k)dk to estimate 
(TrS^) provides a systematic over-estimate; for the lattice we considered, it over-estimates 
the actual magnetic energy by 3.5%. But /e gets the electric energy exactly right. With 
these remarks in mind, we will use fEt to estimate the occupancy and will use /a as a 
secondary estimator, with some caution about the meaning of its infrared behavior. 

3.2 Approach to scaling 

We have performed a number of classical gauge field evolutions on lattices of size up to 
256^ and for times in excess of Qt = 10^. Provided that we concentrate on that part of an 
evolution where apmax < 1, we consistently see the approach to scaling expected from our 
kinetic theory study. To illustrate this. Figure ^ shows g'^Nc f versus p/Q at a number of 
times, first in absolute units and then building in the assumption of scaling as described in 
Eq. ( p.lOP and Eq. ( |2.11| ). We see that, after rescaling / and p to remove their dominant 
time scaling, the occupancy takes the same form quite accurately over a wide range of 
times. 



Qa = 0.1357 



Qa = 0.1357 




Momentum p 



a 



10 w 

1 1 

0.1 

' 0.01 y 

0.001 w 
0.0001 = 



T" 




t/a = 


100, Qt = 


13.7 


t/a = 


1 000, Qt = 


135.7 


• t/a = 


10 000, Qt = 


1357 


• t/a = 


100 000, Qt = 


13570 


t/a = 


250 000, Qt = 


33925 



1 

Momentum p= p/Q {Q t) ' 




Figure 4: Left: occupancy against momentum at a number of times. Right: the same curves, but 
scaling out the dominant time dependence oi p/Q cx {Qt)^ and g'^Nc / oc {Qt)~ . 



How quickly do the gauge fields approach this scaling solution? To study this, we look 
at the same initial conditions used in Figure |l|. Figure ^ shows the occupancies and 
fEt^ starting at an early time where the very different initial conditions are evident and 
continuing until they have approached a common behavior. The spectra agree at the same 
time that the values of Pmax converge in Figure |, of order Qt = 60. 

Looking at Figure ^ and Figure ^, it appears that the scaling solution is not completely 
time independent, but shows some weak time evolution. First, the large-momentum, low- 
occupancy falloff appears to evolve with time. Second, the infrared behavior, below the 
characteristic scale Pmax, also shows some evolution. The evolution at low momentum is 
real and we will look at it in more detail in the next subsection. But the evolution at 
large momentum is a lattice artifact, caused by higher dimension operators which become 
important at momentum scales p ^ 1/a. To illustrate this, we plot the UV falloff at 
several times Qt and lattice spacings apmax in Figure The figure shows clearly that the 
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1283 SU(2) Qa=.42 
Time Qt^8 vs. Fig lA 
Black solid: fg^ 
Blue dotted: f. 
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_ 1283 3u(2) Qa = .42 
Time Qt^64 vs. Fig IB 
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Blue dotted: f. 





1283 3u(2) Qa = .42 
Time Qt^B vs. Fig IB 
Black solid: fj.^ 
Blue dotted: f. 
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Time Qt^24 vs. Fig IB 
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Time Qt^400 vs. Fig IE 
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Figure 5: Approach to the scahng solution for the same 6 initial conditions considered in Figure 
|l|. Initial differences are notable at = 8 but already small at Qt — 24 and gone by Qt = 64. 
Differences in the infrared in the last line are of order statistical error. 



differences between curves are associated with their taking different values of apx\ 
with their being at different times QtJ 
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Figure 6: Ultraviolet falloff from several evolutions, at different times and different lattice spacings. 
The UV behavior is mostly sensitive to apmax; at equal apmax it is almost Qt independent. 



^We should clarify that the momentum in the figure is really the square root of the free lattice dispersion 



which for our Wilson lattice action is P = ^ "Y^^ so? i^^^^- 
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3.3 Evolution of the scaling solution 

Now we examine in more detail whether there is time evolution in the infrared part of the 
scaling solution. We saw in Figure ^ that distinct initial conditions approach a common 
scaling solution by Qt = 60. Our focus will be on time scales this long or longer. To 
do so we performed several independent evolutions on large (256^) and relatively coarse 
(Qa = 0.422) lattices, averaging and obtaining ensemble error bars. For times longer than 
Qt = 500 we also performed simulations at Qa = 0.298. 

The results are displayed in Figure The figure first compares fEt (black solid lines) 

—4 

to /a (blue dashed lines), in the top two figures. To test for a p~3' infrared scaling behavior, 

4 4 ~ 

we have multiplied the occupancies by ps and plotted f{p) in the figure. 

-4 

The lower right plot focuses on It shows that /a scales as /a oc p 3 from about 
p = 0.4 out to the infrared scale where the occupancy hits the value / = 6/g'^Nc, indicated 
by the red dashed lines. Then it saturates at this value. The value Q/g^Nc is the same 
as the maximum value we found in the equilibrium study, so this appears to indicate that 
the occupancy has become nonperturbatively large and encountered magnetic screening. 
We can estimate the total particle number in nonperturbative infrared magnetic fields by 
integrating fA{k)d?k starting where fAik) > 'i/g'^Nc- We find that the fraction of particle 
number in the condensate, 

ncondcnsatc ^ IfA>4^/ig''N,) ^"^ fAdk 

ntotal J k'^fAdk 

is 4.3%, 1.9%, 1.1%, 0.75%, and 0.54% for Qt = 32, 72, 156, 299, and 544 respectively. 

Now consider the electric occupancies. The lower left figure compares fEt with Je- It 
also presents an estimate of the Debye screening scale, obtained by evaluating 

the quantity which would perturbatively give m^. This determination of is indicated 
for each evolution by a (magenta) vertical bar, with the rightmost (large-fc) bar representing 
the earliest time and the leftmost (small-A;) bar the latest time. As we have seen, ttid is 
the scale where we would expect Je and fEt to start diflFering from each other significantly, 
and this is indeed the case. Since fE is established assuming there are only two occupied 
polarizations, while below there are really three, it is fEt which should be used as an 
estimate of occupancy. The fact that fE rises above fEt could lead to an incorrect inference 

-4 

that at least the electric occupancies rise more steeply than p s in the infrared. In fact 
fEt rises in the IR for Qt = 544, but not as strongly as we would find using fE- 

The use of fEt to estimate occupancies is also not really correct below the scale Wpi = 
"i-d/\/3) because Eq. (J^^) and Eq. ( |3.4] ) assume free dispersion of plasmons, while the 
plasmon dispersion is strongly modified in this region; the factor l/\p\ in Eq. dJ) should 



presumably be replaced by 1/yp^ + Wpj. (This is also not quite right because part of the 
electric field strength resides in the Landau cut. We have no good way of separating the 
pole and cut contributions, but for k < ujpi one expects the pole (plasmon) contribution to 
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Figure 7: Occupancies at intermediate times, explicitly scaling out the dominant time dependence 
and the p'S' momentum dependence, based on 7 independent evolutions on 256^ lattices. Top 
left: electric fst (solid black) versus magnetic Ja (blue dashed) occupancies; the earliest time is 
the bottom curve ai p = 0.03 but the top Ja curve at p = 0.4. Top right: zoom-in of the top-left 
figure. Bottom left: comparison of /^ij and Je (electric occupancies contaminated with longitudinal 
modes, shown in dotted green). The (magenta) vertical bars are estimates of the Debye mass, with 
the earliest time at the right. Bottom right: magnetic occupancies, compared to the value where 
/ — the maximum value observed in equilibrium; the lowest dashed curve is for the earliest 

time, the highest for the latest time. 



dominate, so we will ignore this complication.) If we bravely redefine fEt as 

/ijjplasmon] = ^T'^ob f e^P-^(^?(x)4(0))eoui (3.12) 



and re-plot just /^jjplasmon] against p, we get the occupancies shown in Figure |8|. (The 
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figure adds two lines at later times Qt = 976, 1680 obtained on a finer lattice with Qa = 
.298.) 

fj.^ including o. 



T 1 — I — I — r-r-| r 




P = (p/Q) (Qt)-i/^ 

Figure 8: Occupancies fst assuming dispersion corrected by the plasma frequency. Data is the 
same as in Figure ^ but with two later times (at Qa = .298) added. 

What Figure ^ shows is that the apparent rise in , found in Figure 0, arose because 
the occupancies were being mis-evaluated, not taking into account the dispersion of the 
plasmons which /^^ is describing. When one takes this into account, one finds for Qt > 200 

-4 

that there is a clean scaling window with p s scaling behavior, and then a falloff at momenta 
below the plasma frequency. 

In summary, as Qt increases, a scaling window opens between the scale Pmax and 
an infrared scale where occupancies become nonperturbatively large and the physics of 

—4 

plasmons is important. This scaling window shows / oc p s . After correctly identifying 

-3 

the relation between electric field correlators and occupancy, we find no evidence for p 2 
scaling at any time. 

4. Are there condensates? 

Recently Blaizot et al have argued that the evolution of classical Yang-Mills theory, of 
precisely the sort considered here, may lead to the formation of condensates of gluons in 
the deep infrared Q. If by condensate one means that the occupancy reaches the scale 
1/5^ in the infrared, then we definitely see a condensate.^ As discussed above, we find that 

4 

f{p,t)^^{Qt)^(^Y where p = Q{Qt)K (4.1) 
g^Nc \pj 

*Using this definition, there is also a condensate in weakly coupled Yang-Mills theory in equilibrium. 
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which is easily solved for the scale where f{p) ~ l/g^Nc- 

for p~Q(Qt)^. (4.2) 

The occupancy integrated at and below this scale is of order f ~ ?jy {Qt)~ ■ Therefore 
the particle number stored in the condensate decays with time as the —6/7 power of time. 
Relative to the total particle number, which scales as {Qt)~ , the condensate makes up 
a fraction of particle number which scales as {Qt)~ . Using the estimates for magnetic 
particle number with Ja > ^/g^N^ in the last section, we find the data fit this trend 
well, with a prefactor of about 0.45 - that is, if we define as "condensate" any modes 
with occupancy fA > ^/{g'^^c), we find the condensate makes up 0.4:5{Qt)~ of the total 
particle number. 

We do not see evidence for a time-independent or long-lived transient population of 
1/g'^ occupancy modes in excess of the time-scaling estimate above. We also do not see 
evidence for a more robustly-defined condensate, in which the occupancy in a very narrow 
momentum range exceeds l/g'^. We believe that such an / ^ 1/5^ condensate is physically 
possible for plasmons, but happens not to occur; and that it is very difficult for such a 
condensate even to occur in magnetic fields. We will present both of these arguments in 
more detail in this section. 

4.1 Electric condensates: decay rate of plasmons 

One possibility is the development of a condensate of plasmons, which are the low mo- 
mentum extension of the conventional quasiparticles. Is there any evidence for such a 
condensate in the preceding results on occupancies? To answer this question, we think it is 
useful to see what a condensate of plasmons would look like and how it would evolve. We 
can do so by artificially introducing a condensate of plasmons into a classical simulation 
during the cascade towards the ultraviolet, examining it with the tools of the previous 
section. 

With this in mind, we have performed three simulations with Qa = 0.42 on 256^ 
lattices in which we stop the evolution at some point in time, fix to Coulomb gauge, add 
a perfectly spatially uniform electric field, re-enforce Gauss' Law^, and then follow the 
evolution. We performed three runs, introducing the electric field at times Qt = 104, 
Qt = 230, and Qt = 440. 

Figure |^ plots the electric and magnetic energy as a function of time after the intro- 
duction of the IR electric field. Unsurprisingly, the electric field energy jumps up at the 
moment when we add the coherent £'-field. The electric field energy then oscillates between 
an almost equal division with magnetic energy and an excess of electric energy. These are 
plasma oscillations. In a plasma oscillation, energy goes back and forth between the elec- 
tric field and a coherent motion of the charged quasiparticles. Since quasiparticles have 
almost equal electric and magnetic energy, when the plasmon is stored in quasiparticles, 

^Gauss' Law is that DiEi = 0. A uniform electric field in Coulomb gauge satisfies diEi — 0, which is 
not the same; therefore the configuration where a uniform E field has been added will not satisfy Gauss' 
Law. We enforce Gauss' Law by projecting to the Gauss constraint using the algorithm of Ref. [H. 
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Figure 9: Electric and magnetic energy as a function of time, when a uniform electric field is 
artificially introduced. Vertical axis is {E"^) and {B'^) relative to their initial average. At time 
(Qi)juinp, a uniform electric field is added, leading to (plasmon) oscillations between electric and 
magnetic energy. 



{E"^) ~ {B"^). The time between maxima (or between minima) of the electric field energy 
is 7r/wpi. We can estimate Wpi by using = m^/S and getting from Eq. ( 3.11 ), or we 
can read Wpi directly off the plot; the two agree at the 10% level. In particular it is clear 
that ojp\ is smaller at later times. 

The amplitude of the oscillations in E'-field energy decay with time. This indicates 
that the plasmons are either scattered to other momenta, or absorbed by number-changing 
processes. ^'^ The occupancy falls by about a factor of 2 in one plasma oscillation {t = 2t: /uj^\ 
or two {E"^) peaks), with the decay occurring faster for the system where we introduced the 
£^-condensate at an earlier Qt. The time scale for plasmon decay is short compared to the 
system age for all of the cases we considered. Therefore we conclude that a condensate of 
plasmons would be short-lived; at any time, any plasmons must be of recent origin rather 
than quanta left over from the initial conditions. 

We have also investigated the occupancies and Ja when this large plasmon con- 
densate is present. Figure |l^ shows the occupancies for our (Qt)spike = 230 simulation at 
two points during the plasmon oscillations: when the E'-field energy has its first minimum 
and all energy is in quasiparticles (left), and when the £'-field has its next maximum and 
the energy is in an E-^ield condensate (right). The dotted lines in the figure are the oc- 
cupancies at the same time but in a simulation where the S-field spike is not added, for 
comparison. The figure shows that, when {E"^) ~ {B"^), the plasmon energy is stored in 



^"if the plasmons were scattered to momenta differing by <, Upi from the original value, they would 
lose phase coherence with the fc = plasmon, but would continue to have a time-averaged {E^) > (-B^). 
This would prevent the maxima in (B^) from almost reaching the {E"^) minima. But we observe that they 
continue to almost meet; so we believe the plasmon damping is dominated by absorption and scattering to 
momenta k > cjpi. 
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Figure 10: Left: comparison of electric occupancy Je (black) and magnetic occupancy Ja (blue) 
a quarter-oscillation after an i?- field spike is added (first moment when (i?^) ~ (-8^)). Dotted lines 
are occupancies at the same time, for a run without the i?-field spike, for comparison. Right: the 
same but a half-oscillation after the spike, when {E"^) has its next peak. 



a redistribution of the quasiparticles. When {E'^) has its peak, there is a condensate of 
soft i?- fields, with the occupancy reaching 3 orders of magnitude higher than without the 
plasmon condensate. This is a real condensate! By comparison, we can conclude that there 
is not a condensate of plasmons in the normal simulation. 

What this investigation shows is that a condensate of plasmons is possible in principle. 
In particular, an intense and uniform E-field does not create an exponential instability and 
break up in a time scale r ~ l/\/gE Instead, such a condensate oscillates between 
manifesting as an E'-field and as a coherent motion of partons. (Note that the right graph 



in Figure 10 is at a time after the left graph in the figure, and shows that a large -E-field 
occupancy regenerates from the coherent motion of quasiparticles.) 

The evidence for such a condensate in an ordinary simulation would come in two parts; 
an excess of electric over magnetic energy, and a spike in the small-A; occupancy Je- A 
slight excess of {E'^) over {B'^) is observed, but it is already below 2% for Qt = 50 and 
it decays thereafter. Also, such an excess is expected because of ordinary perturbative 
interactions, which should lead to an excess of order m^/p^ ~ {Qt)~ (the scaling we 
observe). We also definitely do not observe an infrared spike in Je, as shown in Figure 
and Figure 0. So while an electric or plasmon condensate is possible, we do not see one. 



4.2 Magnetic condensates: Nielsen-Olesen instabilities 

What about condensates in the magnetic sector? We will argue here that such condensates 
are difficult to define or identify in Coulomb gauge. But we will also show that genuine 
condensates, in the sense of magnetic field occupancies greatly in excess oil/g^ for a narrow 
range of A;-vectors, are highly unstable in a nonabelian theory. 
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First we explain why it is difficult for mag- 
netic occupancies to be highly coherent in the 
sense of occupancies exceeding /a ~ 1/5^, or 
equivalently {AA{k)) ~ k'^/g'^- Suppose that 
there is a gauge field configuration with in- 
frared fields of size A{x) ~ l/{gL), with L 
some length scale. Consider Figure |ll] and ask 
if the field can be coherent between points A, 
B, and C, separated by distances of order L. 
We might say that the field at A and at B 
are "the same" if the parallel transport of the 
field at A, along the path AB, is the same or 
similar to the field at B. Similarly, the field 
is the same between A and C if paral- 
lel transported along path AC is the same as 
A(C). But because the field is large, the Wil- 
son line along AB and then BC will differ by 
an order- 1 group element from the Wilson line 
along AC - roughly, by g times the flux of B field going through the loop, which is order 
1. Therefore, the parallel transport of the field at A, first along path AB and then along 
path BC, cannot agree with the same field parallel transported just along path AC. So a 
comparison of the field at B and at C, using the path BC, will show that they are very 
different. Therefore there is no sense that the field is really coherent on scales as large as 
or larger than L. 

Since there is no gauge invariant sense in which a gauge field of amplitude |A| ~ ^/gL 
can be coherent on scales ^ L, it would be surprising - and probably just a gauge choice 
artifact - for A to be coherent on longer scales. When a field is random over scales longer 
than L, its Fourier transform has roughly uniform power over momenta less than k ~ 1/L. 
If the real-space field strength is ~ l/{g'^L'^) and this is distributed over Fourier modes 
with k ~ 1/L, then the Fourier-space correlator will be {A'^{k)) ~ k/g'^, corresponding to 
an occupancy /^i ~ for /c < 1/L (see Eq. (|3.lD ). 

A loophole in the above reasoning is the case where A points in the same color direction 
throughout a large region of space - that is, if A happens to be approximately abelian (or 
more generally, if A commutes with the field strength through any arbitrary Wilson loop) . 
This is very nongeneric, since A is specified by 3(A'^— 1) numbers (each color can point in 
a different space direction), but a gauge transformation contains only (A'^^— 1) independent 
elements. Therefore 2{N^—1) elements of A must coincidentally vanish for it to be nearly 
abelian. Nevertheless it is useful to consider what happens in this special case, where 
there certainly can be coherence on large scales. So consider the case where only the 
component of A is nonzero, and B is uniform in space. It is quite easy to consider this 
possibility on the lattice, and we have done so, adding roundoff error sized fiuctuations in 
other field modes. We observe the exponential growth of electric field fluctuations; we plot 
{B"^) and {E"^) as a function of Qt in Figure and /a and Je as a function of p/Q at a 




Figure 11: Illustration of the obstruction to 
a magnetic condensate: if parallel transport 
of the field at A, along path AB, agrees with 
the field at B; and transportation along path 
AC agrees with the field at C; then since the 
Wilson loop is very nontrivial, transportation 
from B to C must pick up a large color rota- 
tion. 
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series of times in Figure |T^. (In both figures we set the zero of time to the first moment 
when {E"^) = {B"^).) Figure |T2| shows that the electric field energy grows exponentially until 
it actually exceeds the magnetic energy; then there are brief, highly damped oscillations, 
ending in an almost equal energy partition between E and B. The lefthand plot of Figure 
13 shows occupancies during the process of electric field growth, with successive curves 
at successive times separated by about Qt = 0.8. The righthand plot shows the process 
of field saturation at a series of times around t = 0. The figure shows that the magnetic 
condensate (the spike at the earliest time) disappears; then there is some chaotic dynamics 
with a very short-lived high IR occupancy of electric fields. But these fields are mostly 
gone by Qt = 3.2, when the configuration already resembles the occupancies during the 
cascade. 
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Figure 12: Electric and magnetic energy fractions 2{E'^)/{{E^) + (B^)), 2{B^)/{{E^) + (B^)) as 
a function of time during Nielsen-Olesen instabilities. Left: linear plot. Right: log-linear plot. 



The reason for the exponential growth of electric fields in the presence of a uniform B 



field is the Nielsen-Olesen instability. As shown by Nielsen and Olesen in a uniform 
B field there are modes whose amplitude will grow at the rate^^ 



oc exp(-|-7t) , 



7 



gB. 



(4.3) 



The occupancy and energy go as E'^ and grow at twice this rate. This is precisely the growth 
rate we observe. The arguments of Nielsen and Olesen apply whenever the magnetic field is 
coherent over a distance scale large compared to i?Lamor = 1/ VgB. This essentially forbids 
magnetic fields with such coherence as to have /a ^ l/^^; any such magnetic field would 
be exponentially unstable and rapidly dissolve. 



^^Here is a lightning quick derivation of this result. In the presence of a magnetic field, charges follow 
Lamor orbits, causing transverse momentum to be quantized, fc^ ~ gB{l + 2n), with n = 0,1,2,... 
the Landau level. The energy of an excitation is e = — "igS ■ B + fcy , with —2gS ■ B spin-magnetic 

interaction term. For a spin-1 particle, ~2gS ■ B = ±2gB depending on polarization. The minimum of e is 
obtained when n = 0, fey =0 and —2gS ■ B = —2gB. In this case the energy is e = ^/~gB, corresponding 
to exponential growth with exponent exp{t^gB). 
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Figure 13: Growth of Nielsen-Olesen instabilities for a uniform abelian i3-field. Left: the electric 
occupancies (dashed lines) grow in the presence of the magnetic condensate (solid, spiky line). 
Right: the B-field condensate disappears, there is briefly a very large i?-field occupancy in the 
infrared, and then the configuration becomes the first stage of the usual cascade to the ultraviolet. 



A curious feature of Figure 13 is that, even though we started with a perfectly uniform 
and abelian gauge field configuration, the process of finding Coulomb gauge has interpreted 
it as a nonuniform and nonabelian configuration; the occupancy /a in Figure 13 has its 
peak at A; ~ Q rather than at A; = 0. It is still recognizable as a condensate; using 
our former definition that "condensate" means Ja > ^/{g^^c)-, we find about 90% of the 
particle number is in a condensate. The fact that the field is highly coherent is stored 
in the very peculiar structure of the inferred occupancies, with a sharp spike reaching far 
above / = 6/g'^Nc, the limit we have observed in thermal and stochastic cascading systems. 
We believe that a genuine magnetic condensate would appear as a similar, structured and 
sharp feature in the magnetic occupancy - a feature which does not appear at any point 
in our simulations except those which start with a coherent B-&eld. 

To summarize, we directly observe the Nielsen-Olesen instability; whenever the mag- 
netic field B is coherent on a scale larger than the Lamor radius of charges in the field, 
-RLamor = ^/VdB, there are exponentially unstable modes, growing on a time scale 7 = 
\/ gB, which consume the -B-field and turn into a stochastic collection of excitations. This 
physics should prevent fields with k < \/gB, corresponding to occupancies larger than 
roughly /a Z 2TT/g^. 



5. Discussion 



Classical Yang-Mills theory, initially occupied in the infrared, sees a cascade of energy 
towards the ultraviolet, due to the infinite phase space available there. Simple arguments, 
which suggest Pmax ~ QiQi)^ and g^Nc f ~ {Qt)~ , are borne out, and in fact, scaling out 
this dominant behavior, we find that the occupancy f{p, t) rapidly approaches a scaling 
solution. 
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The infrared dynamics at intermediate times {Qt ~ 100) are complicated because the 
screening and nonperturbative scales are not well separated from the scale Pmax* 

At the 

_2 

scale mD ~ Pma,^{Qt)~ the electric field takes on significant longitudinal occupancy, and 
one should use only its transverse components to estimate the quasiparticle occupancy, 
to avoid miscounting. Below around Wpi = mD/\/3 the dispersion of plasmons should 

_3 

also be considered. Taking these effects into account, we find no evidence for p~ scaling; 

-4 

when any scaling window exists in the infrared, it always shows exponent p 3 . The occu- 
pancy determined from A-field correlators agrees with that from £^-field correlators until 
around the scale ttt-d, where they start to diverge. The magnetic occupancy /a saturates 
at about Ja = Q/{g^Nc), and the total particle number in such nonperturbatively large 
fields ( magnetic condensate) scales with time as "'^"'j^'^""'^"*^" ~ 

It is possible for electric (plasmon) condensates to exist; indeed, we can artificially 
introduce them into a lattice evolution to investigate their behavior. The plasmons can 
carry occupancy large compared to 1/5^, but in practice they do not. And the damping 
rate for plasmons is fast compared to the system's age. 

Genuine magnetic condensates are unstable due to the Nielsen-Olesen instability; we 
verify the growth rate of fluctuations in the presence of a perfectly coherent -B-field, 7 = 
yJgB. The high-occupancy infrared magnetic field observed during the scaling cascade 
to the ultraviolet better resembles the nonperturbative incoherent i?-fields in equilibrium 
weakly-coupled Yang-Mills theory. 

It would be interesting to solve the Boltzmann equations for Yang-Mills theory, pre- 
sented in Reference 0, directly, to compare against the numerical scaling solution found 
here. 
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A. Equilibrium i?-field correlators 

Consider classical Yang- Mills theory with some UV regulator (such as the lattice) . Because 
the theory is UV regulated it has a well defined thermal ensemble, described by the path 
integral 

Z = j V{A,E,) exp ^-^1^1±^^ 6{D,E,) . (A.l) 
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The delta function enforcing Gauss' Law is a functional delta, with one such delta at each 
point. We can rewrite it in terms of a Lagrange multiplier, suggestively named Aq |15]: 

5iD,E,) = J VAo exp (l^^^j . (A.2) 

If we perform the Gaussian integration over the -E-field we obtain 

V{Ao,Ei) exp I — j = VAo exp I — 1 , (A.3) 

the standard bare kinetic term for the field. Under renormalization this mixes with the 
identity, generating the Debye mass 

(A^o)' ^IR {DiAof + mlAl . (A.4) 

But this IR behavior should already hold even if we don't integrate out the i?-field. In the 
infrared, the behavior of Z is 

2iR = / V{Ai,Ei,Ao) exp . (A.5) 



Now we can integrate out the held, 

ZiR = y V{M,E,) exp I ^^3 3) \ ^ (A.6) 

neglect the difference between dj and L'j, and use this to evaluate the transverse and 
longitudinal -E- field correlators: 



hk,{EiE,{k)) = , , (A.7) 



T 

1 + k'^/m'^ 

^Vi,{EiEj{k)) = T . {AI 



Combining with Eq. ( p. 51 ) and Eq. (3.6), Eq. ( |3.7D follows. The sum rule is not exact 
because we have approximated DiEi with diEi; but our numerical results in Figure ^ show 
that it is surprisingly accurate. 
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